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Gene segments from other organisms, such as viruses, are detected as foreign and targeted for silencing by 
RNAi pathways. A deep-sequencing map of the small RNA response to repeated transgenes introduced to 
Caenorhabditis elegans revealed that specific segments are targeted by siRNAs. Silencing of the foreign gene 
segments depends on an antiviral response that involves changes in active and silent chromatin modifications and 
altered levels of antisense siRNAs. Distinct Argonaute proteins target foreign genes for silencing or protection 
against silencing. We used a repeated transgene in a genome-wide screen to identify gene disruptions that enhance 
silencing of foreign genetic elements and identified 69 genes. These genes cluster in four groups based on 
overlapping sets of coexpressed genes, including a group of germline-expressed genes that are likely coregulated by 
the E2F transcription factor. Many of the gene inactivations enhance exogenous RNAi. About half of the 69 genes 
have roles in endogenous RNAi pathways that regulate diverse processes, including silencing of duplicated genes 
and transposons and chromosome segregation. Of these newly identified genes, several are required for siRNA 
biogenesis or stability in the oocyte-specific ERGO-1 pathway, including eri-12, encoding an interactor of the 
RNAi-defective protein RDE-10, and ntl-9/CNOT9, one of several CCR4/NOT complex genes that we identified. 
The conserved ARF-like small GTPase ARL-8 is required specifically for primary siRNA biogenesis or stability in 
the sperm-specific ALG-3/4 endogenous RNAi pathway. 
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Small RNAs (18-32 nucleotides [nt]) regulate gene ex- 
pression, chromosome segregation, DNA elimination, 
and viral defense through post-transcriptional and tran- 
scriptional mechanisms. Exogenous RNAi, the RNAi 
response to administered dsRNA or siRNAs, is a viral 
defense pathway; protein factors mediating exogenous 
RNAi are also required for the silencing of viruses (Lu 
et al. 2005; Schott et al. 2005; Wilkins et al. 2005; Felix 
et al. 2011). Factors required for Caenorhabditis elegans 
exogenous RNAi include many conserved proteins, such 
as the endoribonuclease Dicer (DCR-1) (Grishok et al. 
2001; Ketting et al. 2001), the Argonaute protein RDE-1 
(Tabara et al. 1999), and a dsRNA-binding protein, RDE-4 
(Tabara et al. 2002), as well as several effector Argonaute 
proteins. In addition, C. elegans uses proteins shared with 
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plants and fungi but few other animal species for the 
amplification of the RNAi response — most importantly, 
RNA-dependent RNA polymerases (RdRPs) that amplify 
secondary siRNAs (Sijen et al. 2001). Exogenous RNAi 
triggers not only destruction of the targeted mRNA in the 
cytoplasm but cotranscriptional silencing of the same 
target genes, mediated by the NRDE proteins (including 
an Argonaute protein, NRDE-3, which shuttles siRNAs 
from the cytoplasm into the nucleus), to inhibit RNA 
polymerase II elongation and induce deposition of H3K9me3 
chromatin marks on the genetic locus targeted by com- 
plementary siRNAs (Guang et al. 2008, 2010 ; Burkhart 
et al. 2011; Gu et al. 2012). 

The complex machinery of RNAi has regulatory roles 
apart from immunity against newly introduced foreign 
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genetic elements. Endogenous RNAi pathways that si- 
lence a range of genes resident in genomes have also been 
identified in yeast, plants, nematodes, fruit flies, and 
mice. C. elegans endogenous RNAi pathways can be 
distinguished by the distinct Argonaute proteins and 
the length and 5' nucleotides of the siRNAs as well as 
the gene loci from which the small RNAs are derived. 
In oocytes and embryos, the ERGO- 1 -associating 26G 
siRNAs and NRDE-3 -associating 22G siRNAs silence 
recently duplicated genes (Vasale et al. 2010; Fischer 
et al. 2011); the ALG-3 and ALG-4 26G siRNA pathway 
in sperm is required for sperm morphogenesis (Conine 
et al. 2010). CSR-1 -associated 22G siRNAs direct chroma- 
tin modifications on germline-expressed genes (Clay comb 
et al. 2009). Many 22G siRNAs associate with C. elegans- 
specific Argonautes (the so-called WAGOs); these siRNAs 
are produced from many genes, including transposons, 
and direct the silencing of transposons (Gu et al. 2009). 
Several of the endogenous RNAi pathways act in the 
nucleus. The ERGO-1 endogenous RNAi pathway acts 
upstream of the NRDE-3 Argonaute that shuttles endog- 
enous siRNAs into the nucleus. NRDE-3 in conjunction 
with other NRDE proteins directs cotranscriptional si- 
lencing of siRNA target genes and directs deposition of 
the inheritable silent chromatin mark histone H3K9me3. 
The CSR-1 endogenous RNAi pathway is required for 
chromosome segregation (Claycomb et al. 2009) and is 
also thought to confer protection against silencing (Lee 
etal. 2012). 

Small RNA pathways compete for limiting factors 
that are required in multiple pathways (Ambros et al. 
2003; Lund et al. 2011). A few genes that function as 
negative regulators of the C. elegans exogenous RNAi 
pathway have emerged from noncomprehensive genetic 
screens for enhanced response to exogenous RNAi, 
revealing the ERI-1 nuclease (Kennedy et al. 2004), the 
ERI-2/RRF-3 RdRP (Simmer et al. 2002), ERI-3 (Duchaine 
et al. 2006), an activating mutation in DCR-1, the Dicer 
ortholog (Pavelec et al. 2009), the ERI-5 Tudor pro- 
tein (Duchaine et al. 2006), the ERI-6/7 helicase (Fischer 
et al. 2008), ERI-9 (Pavelec et al. 2009), and the Argonaute 
ERI-8/ERGO-1 (Pavelec et al. 2009; Fischer et al. 2011). 
The activity of these genes normally attenuates the re- 
sponse to ingested or injected dsRNA but is also required 
for certain endogenous RNAi pathways. The concomitant 
increase in exogenous RNAi response and decrease in 
endogenous RNAi response may be due to competition 
for particular limiting factors that are shared between 
multiple distinct small RNA pathways,- for example, 
an Argonaute protein (Yigit et al. 2006). Alternatively, 
factors repressed by the endogenous RNAi pathway could 
encode limiting components of the exogenous RNAi 
pathway. 

The closest homologs of many endogenous RNAi 
pathway factors identified in C. elegans fulfill similar 
functions in higher organisms,- e.g., the endogenous 
siRNA biogenesis machinery in the ERGO-1 pathway 
resembles piRNA biogenesis complexes identified in 
flies and mammals, the ERI-6/7 helicase (Fischer et al. 
2011) may be functionally equivalent to its putative 



orthologs Drosophila Armitage (Saito et al. 2010) and 
mouse MovlOLl (Zheng et al. 2010), and the Piwi-like 
Argonaute ERGO-1 could act like Drosophila Piwi and 
mouse MILI and MIWI2. The endogenous siRNAs pro- 
duced by the ERI-6/7 helicase and the ERGO-1 Argonaute 
resemble piRNAs in terms of 2'-0-methylation of the 
3 '-terminal nucleotide, and their potential to trigger 
siRNA biogenesis in trans is similar to piRNAs (Bagijn 
et al. 2012,- Montgomery et al. 2012). 

Transgene silencing is mediated by many of the same 
factors as those needed for RNAi (Tabara et al. 1999; 
Ketting and Plasterk 2000; Vastenhouw et al. 2003; 
Grishok et al. 2005; Robert et al. 2005), and transgene 
silencing is enhanced in mutants that display enhanced 
exogenous RNAi responses (Simmer et al. 2002), arguing 
that the foreign and repeated genetic elements on trans- 
genes are silenced by RNAi. Highly repetitive transgenes 
also accumulate heterochromatic marks concomitantly 
with silencing by RNAi pathways (She et al. 2009; Towbin 
et al. 2010; Arico et al. 2011), suggesting that they are 
subject to transcriptional gene silencing in addition to 
mRNA degradation. Thus, the silencing or desilencing 
of transgenes is an excellent surrogate assay for defects 
in small RNA-mediated silencing. In addition, because 
transgenes bear many of the foreign marks that are also 
associated with viruses, the RNAi response to trans- 
genes may actually be closer to the biology of RNAi as 
an antiviral defense than the rather artificial feeding or 
injection of dsRNA that is used for exogenous RNAi in 
the laboratory. The continued silencing of transgenes 
long after the initial exposure to foreign DNA in trans- 
genesis points to a capacity to mark these genes as foreign 
and regenerate those marks. 

Here we describe a comprehensive analysis of the 
silencing of a chromosomally integrated, multicopy 
transgene. The silencing response to this transgene in 
wild type is similar to antiviral responses: It is mediated 
by antisense siRNAs and antiviral proteins but also 
involves changes in active and silencing chromatin 
marks and multiple Argonaute proteins exerting oppo- 
site silencing effects. In a whole-genome RNAi screen 
for gene inactivations that cause increased silencing of 
a chromosomally integrated, multicopy transgene, we 
identified dozens of gene inactivations that cause in- 
creased transgene silencing. Analysis of transcriptional 
coexpression shows that these genes cluster into four 
groups that appear to be coregulated. Our secondary 
screens for defects in exogenous RNAi show that many 
of the genes identified are negative regulators of the 
antiviral exogenous RNAi pathway. Thus, inactivation 
of these genes, many of which are conserved in humans, 
could promote immunity against viruses. Our secondary 
screens for defects in endogenous RNAi revealed that 
many of the genes identified mediate multiple gene 
silencing pathways. These gene inactivations may iden- 
tify competing pathways that, when inactivated, in- 
crease the capacity of the exogenous RNAi pathway to 
detect transgenes. As examples of genes acting in en- 
dogenous RNAi, we describe three new genes acting in 
the oocyte-specific ERGO-1 endogenous RNAi path- 
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way. We also demonstrate a role for the conserved GTP- 
binding protein ARL-8 in the biogenesis of primary 
siRNAs in the sperm-specific ALG-3/4 endogenous RNAi 
pathway. 

Results 

A C. elegans siRNA immune response to foreign genes 

Foreign genetic segments can be introduced into 
C. elegans as single-copy integrants or as tandemly 
repeated transgenes that are either integrated into a chro- 
mosome or maintained as a separate unstable minichro- 
mosome. Tandemly repeated transgenes, whether inte- 
grated chromosomally or as minichromosomes, are 
silenced in the germline but not in the somatic cells of 
wild type, whereas single-copy transgenes are not as 
routinely silenced in the germline. Tandemly repeated 
transgenes become silenced in somatic cells as well as 
germline cells in most mutants that cause enhanced 
exogenous RNAi [eh mutants) (Simmer et al. 2002; 
Kennedy et al. 2004; Fischer et al. 2008) and can become 
desilenced in particular RNAi-defective mutants,- for 
example, mut-16 or mut-7 (Simmer et al. 2002). How 
RNAi pathways and silencing of repetitive transgenes 
intersect mechanistically is unknown. We analyzed the 
C. elegans response to two multicopy (>100 copies) trans- 
genes: One transgene, mgls30, contains three genetic 
elements from foreign as well as C. elegans sources 
integrated as a multicopy tandem array into a C. elegans 
chromosome — the foreign Aequorea victoria gfp and 
Escherichia coli lacZ genes and an endogenous collagen 
gene (rol-6) (Fig. 1A). These genes are expressed in 
several somatic tissues of wild type: the hypodermis, 
neurons, and the somatic tissues of the gonad. A second 
multicopy integrated transgene that was tested consists 
of gfp under the control of a ubiquitous sur-5 promoter. 

The genetic markers or GFP from these transgenes are 
expressed in the somatic cells of wild type but are 
strongly silenced in the somatic cells of ^ri-6/7/helicase 
or ergo-l/Argonaute-enhanced RNAi mutants. For exam- 
ple, wild type carrying mgls30 shows a rolling phenotype 
conferred by the dominant mutation in rol-6(sul006). 
However, loss of eri-6/7/helicase or the ergo- 1/Argonaute 
causes complete silencing of the mutant rol-6 collagen 
expression from the transgene,- the mRNA encoding the 
mutant collagen is not transcribed or translated, and no 
rolling phenotype is observed (Table 1). 

Deep sequencing of the small RNA repertoire produced 
from the genome of a multicopy mgls30 transgenic strain 
revealed that siRNAs are generated targeting distinct 
segments of the transgene: the promoter, the ORF, the 
3' untranslated region (UTR), and downstream vector 
sequences (Fig. IB). These siRNAs harbor the hallmarks 
of secondary siRNAs produced by RdRPs in response to 
a primary trigger siRNA or piRNA; these are 22 nt in 
length, and the 5' nucleotide is a guanine (Supplemental 
Fig. 1A). 

Paradoxically, even though the rol-6 genetic element is 
far more potently silenced in an eri-6/7 or ergo-1 mutant 



(Table 1), deep sequencing of three silenced transgenes 
[rol-6 ::lin- 14, col- 10 r.lacZ :: lin-14, and sur-5:: gfp) in the 
eri-6/7 mutant (Supplemental Fig. 1B,C) showed a 10-fold 
reduction in antisense transgene siRNAs targeting pro- 
moter, coding, and 3' UTR sequences, including rol-6 
(Fig. 1C; Supplemental Fig. IB), an observation confirmed 
by quantitative RT-PCR (qRT-PCR) (Fig. ID). As de- 
scribed in more detail below, our analysis favors a shift 
in the balance between the loading of transgene siRNAs 
onto a silencing Argonaute versus an anti-silencing 
Argonaute upon loss of ERI-6/7, accompanied by a histone 
methylation-based silencing that is not associated with 
production of siRNAs in the maintenance of the silenced 
state. Sense siRNAs, present at much lower levels than 
antisense siRNAs, are not reduced in the ergo-1 or eri-6/7 
mutants (Supplemental Fig. ID). siRNAs corresponding 
to lim-6::gfp, a transgene expressed mainly in neurons 
that is not silenced in eri-6/7 mutants, are not reduced in 
eri-6/7 mutants. 

Foreign genes trigger an antiviral silencing response 
involving opposing Argonaute proteins and changes 
in chromatin modifications 

To understand the biogenesis and activity of siRNAs that 
target the foreign and C. elegans gene segments on the 
transgene, we assayed known RNAi factors for their 
requirement for transgene silencing in eri-6/7 mutants 
and, conversely, transgene expression in a wild-type 
background. The ERGO-1 /Argonaute acts in the same 
endogenous RNAi pathway as the ERI-6/7/helicase. Loss 
of ergo-1 or eri-6/7 silences the transgene, as measured by 
the expression of the rol-6 genetic marker (the rolling 
phenotype). The exogenous RNAi factor DRH-l/helicase 
is required for transgene silencing in the eri-6/7 mutant 
(Table 1 ), indicating that a transgene dsRNA intermediate 
may be a substrate of the exogenous RNAi machinery 
DCR-l/Dicer, RDE-4/dsRBD, DRH-l/helicase, andRDE-1/ 
Argonaute. These factors are also required for silencing of 
RNA viruses, suggesting that genetic screens for mutations 
that disrupt the silencing of high-copy somatic transgenes 
reveal antiviral defense pathways. 

Of the siRNAs targeting the rol-6 collagen: :lin- 14 3' 
UTR fusion transgene, most target the lin-14 3' UTR. The 
lin-14 3' UTR is targeted by the microRNAs (miRNAs) 
lin-4 (seven sites) and several let-7 family members. We 
tested and rejected the tantalizing possibility that the 
abundance of siRNAs targeting the lin-14 3' UTR is 
triggered by miRNAs. First, no lin-14 siRNAs are present 
in wild- type worms, where the lin-4 miRNA targets the 
lin-14 3' UTR at late larval and adult stages. Second, 
when we aligned all siRNAs to the lin-14 3' UTR 
sequence in the mgls30 transgenic strain, we found that 
siRNAs map over the entire 3' UTR except 3' of the 
poly(A) signal (Supplemental Fig. IE), although of the six 
most abundant lin-14 siRNAs, five overlap or are adjacent 
to sequences targeted by the miRNA lin-4, let-7, or mir- 
241. Third, the microRNA lin-4 and the miRNA-binding 
Argonaute ALG-1 are not required for transgene silenc- 
ing in the eri-6/7 mutant (Table 1). In fact, loss of the 
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Figure 1. Silencing of multicopy transgenes involves small RNAs, multiple small RNA pathways, and changes in H3K9 and H3K4 
trimethylation. (A) Three plasmids constitute the high-copy mgls30 transgene. [B] The mgls30 transgene produces sense (black) and 
antisense (red) siRNAs, as analyzed by small RNA deep sequencing. Indicated are the small RNAs reads in transgenic C. elegans and 
wild- type C. elegans, corresponding to the rol-6 plasmid that is part of the mgls30 transgene. (C) Transgene siRNAs are reduced in eri-6/ 
7/helicase mutants. (D) lin-14 3' UTR siRNAs are dependent on eri-6 /7helicase and nrde-3 /Argonaute. siR-1 is an eri-6/7-dependent 
endogenous siRNA used as a control. (E) Loss of eri-6/7 results in an increase in the H3K9me3 modification on transgene sequences. 
Loss of the HMT set-25 results in a loss of H3K9me3. E01G4.5 is an endogenous ergo-1 and eri-6/7 target. {F) Loss of eri-6/7 results in 
a decrease in the H3K4me3 modification on transgene sequences. (G) Multiple Argonautes and small RNA pathways act on multicopy 
transgenes. [H) More than one-fifth of endogenous ERGO-1-ERI-6/7-NRDE-3 targets are also targeted by the Argonaute CSR-1, 
including E01G4.5. One-fourth of ALG-3/4 targets are also targets of CSR-1. (/) Model of transgene silencing. 
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Table 1. Gene inactivations that cause a change in the somatic silencing of transgenes in wild type or the eri-6/7 mutant 



Gene A 
Inactivation 


Gene A Description 


Gene B tested for 
suppression 


Gene B Description 


Transgene 
expression 


- (wild type) 








expressed 


Endogenous RNAi factors 


ergo-1(tm1860) 


Argonaute 






silenced 


eri-6/7(mg411) 


helicase 






silenced 


eri-6/7 (mg379) 


helicase 






silenced 


eri-6/7(mg379) 


helicase 


glp-l(RNAi) 


no germline 


silenced 


nrde-3(tm1116) 


nuclear Argonaute 






partially silenced 


eri-6/7 (mg379) 


helicase 


nrde-3(tm1116) 


nuclear Argonaute 


expressed 


wago-1(ok1074) 


Argonaute 






expressed 


csr-l(RNAi) 


Argonaute 






silenced 


Exogenous RNAi factors 


rde-1(ne300) 


Argonaute 






expressed; sick 


rde-4(ne299) 


dsRNA binding 






expressed; sick 


drh-l(RNAi) 


helicase 






expressed 


eri-6/7(mg411) 


helicase 


drh-l(RNAi) 


helicase 


partially silenced 


piRNA factors 


prg-1(tm872) 


Argonaute 






partially silenced 


eri-6/7 (mg379) 


helicase 


prg-1(tm872) 


Argonaute 


silenced 


Chromatin modification factors 


set-25(n5021) 


H3K9 methyltransferase 






expressed 


eri-6/7(RNAi) 


helicase 


set-25(n5021) 


H3K9 methyltransferase 


partially silenced 


lin-61(RNAi) 


H3K9Me2/3 binding 






expressed 


eri-6/7(mg411) 


helicase 


lin-61(RNAi) 


H3K9Me2/3 binding 


expressed 


dpy-30(RNAi) 


Set1/MLL complex 






silenced 


set-2(RNAi) 


H3K4 methyltransferase, Set1 






partially silenced 


eri-6/7 (mg379) 


helicase 


dpy-30(RNAi) 


Set1/MLL complex 


silenced 


eri-6/7 (mg379) 


helicase 


set-2(RNAi) 


H3K4 methyltransferase, Set1 


silenced 


MicroRNA factors 


Iin-4(ma161) 


microRNA targeting lin-14 






expressed 


eri-6/7(RNAi) 


helicase 


Iin-4(ma161) 


microRNA targeting lin-14 


silenced 


alg-1(gk214) 


Argonaute 






silenced 


eri-6/7(RNAi) 


helicase 


alg-1(gk214) 


Argonaute 


silenced 


rol-6(su1006) 


rol-6 not on mgls30 transgene 


alg-l(RNAi) 


Argonaute 


rol-6 expressed 



All mutations and RNAi knockdowns were analyzed in mgls30 transgenic animals and scored for a rolling phenotype. (Red) Gene 
inactivations that cause silencing of the transgene; (blue) gene inactivations that suppress silencing of the transgene. 



Argonaute ALG-1 in a wild-type background silences the 
transgene, as does loss of the miRNA factors VIG-1 and 
AIN-1, the GW182 ortholog (Table 1 ; see below). 

Transgene silencing of single-copy foreign sequences in 
the germline, but not of high-copy transgene sequences in 
the germline, involves piRNAs that target foreign se- 
quences and thus requires the PIWI protein PRG-1. High- 
copy transgene silencing in the somatic cells of the eri-6/ 
7 mutant does not require the PIWI protein PRG-1 (Table 
1). In addition, a functional germline is not required for 
silencing in the eri-6/7 mutant: Transgenes are silenced 
in sterile glp-l(RNAi) animals (data not shown). These 
data indicate that the silencing response to somatic high- 
copy transgenes differs from the response to germline- 
expressed low-copy transgenes,- germline low-copy trans- 
genes are silenced through a response similar to the 
silencing of transposable elements, involving PIWI pro- 
teins and piRNAs. 

The silencing of the rol-6 transgene is suppressed by 
loss of NRDE-3, an Argonaute acting in nuclear RNAi 



(Table 1). This suggests that transgene silencing may 
proceed through a nuclear silencing pathway. NRDE-3 
acts in a nuclear RNAi pathway upstream of the inactive 
chromatin H3K9 trimethylation mark in somatic tissues 
(Burkhart et al. 2011). We tested whether histone H3 
methyltransferases are required for transgene silencing in 
the eri-6/7 mutant. Whereas two members of the Setl/ 
MLL complex (SET-2 and DPY-30) that deposits active 
H3K4me3 marks are not required for silencing, SET-25, 
a histone H3K9 methyltransferase, and LIN-61, an MBT 
protein that preferentially binds H3K9me2 and H3K9me3, 
are required for silencing in somatic cells of the mgls30 
transgene in the eri-6/7 helicase mutant (Table 1). Gene 
inactivation of set-2 or dpy-30 in wild type causes 
silencing of the mgls30 transgene, suggesting that the 
active chromatin marks conferred by these histone meth- 
yltransferases are essential for expression of foreign 
genetic elements from tandem arrays (Table 1). The 
finding that the NRDE-3/Argonaute and SET-25/H3K9- 
HMT are required for transgene silencing in the eri-6/7 
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mutant suggests that transgene chromatin may be modi- 
fied to silence transgene expression. We performed chro- 
matin immunoprecipitation (ChIP) experiments with a 
wild- type mgls30 transgenic strain, an mgls30 transgenic 
eri-6/7 strain (silenced transgene), and an mgls30 trans- 
genic set-25 strain (partially silenced transgene) using 
antibodies for H3K9me3 and H3K4me3. H3K9me3 is 
found at inactive promoters and gene bodies,- in C. elegans, 
it is also found on active genes present in heterochroma- 
tin. H3K4me3 is found on active and inactive promoters, 
including those of genes that may become active, and also 
on the first exon of active genes. 

The silent mark H3K9me3 is increased on transgene 
sequences in the eri-6/7 mutant and reduced in the set- 
25/HMT mutant (Fig. IE; Supplemental Fig. IF). The 
E01G4.5 gene is a known target gene of the ERGO-1- 
ERI-6/7 endogenous RNAi pathway that is desilenced in 
the eri-6/7 mutant; it shows a reduction of H3K9me3 in 
the eri-6/7 and set-25 mutants. 

The H3K4me3 modification is found on the rol-6/ 
collagen exon in a wild-type transgenic strain (Fig. IF; 
Supplemental Fig. 1G) ; in the eri-6/7 mutant, this mod- 
ification is reduced. Similar results were found for the 
E01G4.5 gene. 

These data show that loss of ERI-6/7 leads to changes in 
both active and inactive chromatin marks. As described 
below, we found that loss of the anti-silencing Argonaute 
CSR-1 results in transgene silencing (Table 1; Fig. 1G). We 
hypothesize that transgene siRNAs are bound by both the 
silencing Argonaute NRDE-3 (directing H3K9 trimethyla- 
tion by SET-25) and the anti-silencing Argonaute CSR-1 
(possibly directing H3K4 trimethylation by SET-2/DPY-30). 
CSR-1, ERGO-1, and NRDE-3 also act in endogenous 
RNAi. In an analysis of endogenous RNAi target genes, 
we compared CSR-1 target genes with ERGO-1 target 
genes and found that for more than a fifth of ERGO-1 
target genes, the siRNAs are bound by not only the 
silencing NRDE-3/Argonaute but also the anti-silencing 
CSR-1 /Argonaute (Fig. 1H). This is also true for the gene 
E01G4.5, a gene with two paralogs in the C. elegans Bristol 
N2 genome, all of which are absent in the genome of another 
C. elegans isolate, CB4856 from Hawaii (Maydan et al. 
2007), but present in high copy in other Caenorhabdites: 
Caenorhabditis remanei and Caenorhabditis brenneri. 
E01G4.5 and its paralogs encode proteins with a viral 
DNA polymerase B2 domain, suggesting that they may 
be integrated DNA viruses that continue to be detected 
and repressed by the ERGO-1 -ERI-6/7 newly acquired 
foreign gene surveillance pathway. 

This analysis of multicopy transgene silencing shows 
that these transgenes produce dsRNAs that are substrates 
of the exogenous RNAi and antiviral defense pathway 
mediated by the DRH-1 helicase. ERI-6/7 is required for 
the efficient production or stability of 22G transgene 
siRNAs, although these ERI-6/7-dependent siRNAs are 
not required for silencing, as loss of ERI-6/7 results in 
increased transgene or exogenous RNAi silencing, pre- 
sumably through the release of a limiting factor acting in 
the downstream steps of the silencing pathway. Trans- 
gene siRNAs may also be routed into the CSR-1 endog- 



enous RNAi pathway, a pathway that is thought to confer 
protection against silencing to its target genes (Fig. II; Lee 
et al. 2012). The balance between siRNA loading onto 
NRDE-3 versus CSR-1 may shift upon loss of ERI-6/7. 
Downstream from transgene siRNA generation, the nu- 
clear Argonaute NRDE-3 is strictly required for silencing 
and directs H3K9 trimethylation by the his tone H3K9 
histone methyltransf erase SET-25 and the H3K9me3- 
binding MBT domain protein LIN-61, a protein that 
interacts with SET-25 (X Wu and G Ruvkun, unpubl.). 
This is accompanied by a loss of H3K4me3 marks that 
may have been deposited by SET-2/DPY-30, possibly 
directed by the Argonaute CSR-1. 

Identification of gene inactivations that enhance 
the silencing of foreign genetic elements 

To identify additional genetic components that mediate 
negative regulation or pathways that compete with the 
silencing of repetitive foreign DNA and possibly viruses, 
we performed a genome-wide RNAi screen for gene 
inactivations that enhance the silencing of the rol-6:: 
lin-14 fusion gene on the same mgls30 transgene charac- 
terized above. The mgls30 transgenic strain shows a roll- 
ing phenotype caused by expression of the rol-6(sul006) 
mutant transgene (Supplemental Fig. 2A). The lin-14 3' 
UTR is a target of the lin-4 miRNA and the let- 7 family of 
miRNAs, but normal miRNA activity of wild type is not 
sufficient to silence the multiple copies of the rol-6 ::lin- 
14 3' UTR fusion genes on the tandemly repeated trans- 
gene array. The RNAi screen for a nonrolling phenotype 
was expected to reveal gene inactivations that increase 
the intensity of RNAi, as exemplified by eri mutants, or 
increase the silencing by the miRNAs that target the lin- 
14 3' UTR present on the transgene. 

Transgenic animals were exposed to dsRNAs for two 
generations and scored in the second generation for 
suppression of the rolling phenotype of the parent strain 
(Supplemental Fig. 2B). Ninety percent of the genes in the 
genome that are present in the combination of feeding 
RNAi libraries constructed by Ahringer and Vidal (Boulton 
et al. 2002; Kamath et al. 2003) were tested in duplicate in 
the primary screen. Gene inactivations that repressed the 
rolling phenotype of the mgls30 transgene were retested in 
quadruplicate. Positive hits were tested for genetic in- 
teractions with the rol-6 collagen gene or another cuticular 
collagen gene, sqt-1, using chromosomally localized mu- 
tations,- gene inactivations that suppressed rol-6 and sqt-1 
mutant phenotypes and thus correspond to suppressors of 
these collagen mutations were excluded from further 
analysis. The 69 remaining gene inactivations cause a sup- 
pression of the rol-6 rolling phenotype in >35% of the 
animals in multiple retests (Table 2; Supplemental Fig. 2C) 
only when rol-6(sul006) was expressed from a transgene. 
More than two -thirds of the genes identified in this screen 
are conserved in humans (Fig. 2 A,- Table 2). In addition, 
whereas genetic screens for RNAi pathway components 
have mostly identified gene inactivations that do not 
cause developmental defects, such as arrest of develop- 
ment or sterility, 60% of the genes identified in this RNAi 



GENES & DEVELOPMENT 2683 



Fischer et al. 



Table 2. Sixty-nine gene inactivations that cause enhanced silencing of the mgls30 transgene 



Gene 


Percent silencing 


Description 


Human ortholog 


adr-1 


98% 


dsRNA-specific adenosine deaminase 


AD ARB 2 


U^UUo. 1 


yo /o 


dsRNA-specific adenosine deaminase 




sl-mr OO 


QQ O/ 

yo /o 


Adaptin; clathrin adaptor complex, |x2 subunit of AP2 complex 


A TDO A/T 1 


ergo-1 


98% 


Argonaute protein 




eh- 7 


no o/ 
98 To 


Superfamily I helicase 


IGHMBPZ/ MU V 1 U 


inp-5 


QQ O/ 

yo /o 


Splicing factor exon— exon junction complex protein 


DTVTDC1 


lin-o4 


noo/ 
9o /o 


DNA-binding protein; recruits DRM complex 


T TTVTC A 

LiJNd4 


im-y 


yo /o 


Retinoblastoma pathway protein 


T TTVTO 


ide-4 


yo /o 


dsRNA-binding protein; interacts with Dicer 




eri-6 


93% 






csi-1 


91 % 


Argonaute protein 




aat-z 


QQ O/ 

oy /o 


Secreted metalloprotease 


AlJAlVl 10/ 




QQO/ 

oy /o 


RNA helicase 


UUAoA 


adi-2 


OQO/ 

OO /o 


dsRNA-specific adenosine deaminase 




isw-1 


OO /o 


SWI/SNF-related regulator of chromatin 


C AA ADfA 1 


eii-12 


QCO/ 

OD /o 






Art-rT on 


OJ /o 


BED Zn ringer protein 




ell-1 


OD /o 


RNA polymerase II elongation factor 




ntz-i 


OD /o 


nzA.Z: nistone nz variant 


nzArZ, 


TO A A Q A 


7QO/ 

/y /o 


Zn finger 




stc-1 


77°/ 
/ / /o 


ATPase similar to HSP70 superfamily; ER chaperone 


nolAlo 


suro-1 


77% 


Carboxypeptidase 


CPA2 


ntl-y 


"7/" O/ 

/OTo 


IN (J I -like, CCR4/NOI complex component 


CNUI9 


attj-o 


73 o/ 
/o /o 


AT hook transcription factor 




eKi-i 


/U /o 


Tudor domain protein 




usp-5 


7no/ 
/U /o 


Ubiquitin-specific protease 


T TC "DC 


n/-7V 0 
pCLl-d. 


/:qo/ 
oy /o 


Protein disulfide isomerase 


D/ITTD 

r4rli5 


/no C A IZ O 


68% 


Nuclear export receptor; exportin 


XP07 


prmt-1 


/no/ 
OO To 


Protein arginine methyltransferase 


hLMR 1 1 L2 


COD A ~7 / 

rZyA/.o 


66 /o 


M-phase phosphoprotein 6; exosome associated, RNA binding 


MrnUbrno/MJVLro 


n A Z? 7 
loyKjlUAK. / 


DO /o 






lin-35 


65% 


• 

Retinoblastoma protein (Rb) 


RBL2 


jm-Z 


64 To 


Actin-binding cytoskeleton protein, filamin 


FLNC 


mes-2 


6Z /o 


SET domain protein; Polycomb-like repressive complex PRC2 




sdha-2 


62% 


Succinate dehydrogenase, flavoprotein subunit 


SDHA 


sams-4 or sams—3 


6U /o 


S-adenosylmethionine synthetase 


IVlAl ZA 


ain-o 


coo/ 
Do /o 


DEAD-box helicase 


Irlrll 


tcei-1 


rno/ 

Do /o 


Transcription elongation regulator 




COl-yu 


C/£ 0/ 
D6 /o 


Collagen 


r^fAT O Q A 1 

LULZoAl 


ent-1 


C /l o/ 

54 /o 


Equilibrative nucleoside transporter 


CT POfl A Q /"U "C TVTTQ 

oLCZyAo/nEJN 16 


arl-8 


54% 


Arf-like small GTPase 


ARL8B 


R186.7 


53% 


Pinin 


PNN 


mel-32 


/IQO/ 

4y /o 


Glycine/serine hydroxymethyltransf erase 




T1 1 EO 7 


/IQO/ 

4y /o 






aps-2 


/IQO/ 

4o /o 


Adaptm,- clatnrm adaptor complex, cr2 subunit of AP2 complex 


A DO C 1 

ArZo 1 


dpy-30 


48% 


Histone H3 (Lys4) methyltransferase complex (Setlp) 


DFYdU 


ECV/^O "7 


/IQO/ 

4o /o 






exos-1 


^o /o 


Exosome component 


IjAUjV^I 


TOAf^A O 


/I C o/ 
4D /o 






r~<AQ TJA A 

L l 4orio.4 


/ico/ 
4D /o 






r3bb8.4 


44% 


Metridin-like ShK toxin 




cogc-5 


43% 


Conserved oligomeric golgi complex subunit 


COG5 


ain-1 


43% 


Homolog of TNRC6A/GW182; AGO interactor 


TNRC6A 


cpar-1 


43% 


CENP-A inner kinetochore histone H3 variant 


CENPA 


T03F6.4 


43% 






Y48A6B.7 


43% 


Cytosine deaminase 


ADAT2 


npp-12 


43% 


Nuclear pore complex, gp210 component 


NUP210 


ubc-20 


42% 


Ubiquitin-conjugating enzyme E2 K 


UBE2K 


mes-3 


41% 


Member of Polycomb-like repressive complex PRC2 




cutl-20 


39% 


Cuticulin 





(continued) 
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Table 2. continued 



Gene 


Percent silencing 


Description 


Human ortholog 


vig-1 


39% 


RNA-binding protein RISC complex, homolog of Drosphila VIG 


HABP4 


F37A4.1 


39% 


Abhydrolase 


ABHD16A 


hgo-1 


38% 


Homogentisate 1,2-dioxygenase 


HGD 


puf-8 


38% 


PUF (Pumilio/FBF) domain- containing 


PUM2 


C41D11.4 


38% 






apa-2 


36% 


Adaptin; clathrin adaptor complex, a subunit of AP2 complex 


AP2A1 


C25F9.11 


36% 


G-protein-coupled receptor 


GPR139 


ccr-4 


36% 


CCR4, CCR4/NOT complex component 


CNOT6L 


nurf-1 


36% 


NURF chromatin remodeling complex, subunit NURF301/BPTF 


BPTF 



Inactivation of these genes results in transgene silencing, scored as the percentage of nonrolling animals (percent silencing). Human 
orthologs are indicated. Orthologs in bold have been identified in a comprehensive analysis of the human mRNA interactome (Castello 
et al. 2012). 



screen are essential for viability (Fig. 2A). These genes can 
be detected in functional genomic screens such as these 
that assay for transgene silencing in sterile or arrested 
animals and/or use a shorter gene inactivation period 
than we used for these particular gene inactivations. 

The identification of several gene inactivations known 
to cause increased transgene silencing and/or enhanced 
RNAi, such as the helicase gene eri-6/7 (Fischer et al. 
2008), the Argonaute gene ergo-1 (Fischer et al. 201 1), the 
ADAR genes adi-1 and adi-2 (Knight and Bass 2002), and 
the DRM complex component lin-35/Rb (Kelly and Fire 
1998; Hsieh et al. 1999), validates the ability of the screen 
to discover the previously identified eri genes and regu- 
lators of expression of the RNAi machinery (Wu et al. 
2012). The complete set of genes identified by this screen 
was analyzed for known and predicted interactions as 
well as transcriptional coexpression. We identified mul- 
tiple components of several known protein complexes 
that interact with chromatin, DNA, or RNA, including 
DRM complex genes (lin-9, lin-35, and lin-54), CCR4/ 
NOT complex genes (ccr-4 and ntl-9), PRC2 complex 
genes (mes-2 and mes-3), miRISC complex genes [alg-1, 
vig-1, and ain-1), and NURF complex genes (nurf-1 and 
isw-1). Analysis of interactions between the 69 identified 
genes and proteins yielded a compact functional cou- 
pling subnetwork (Alexeyenko and Sonnhammer 2009) 
with more links (203 links) than the set of 69 previously 
characterized known C. elegans small RNA pathway 
genes (108 links) (Supplemental Fig. 2D). The newly 
identified genes are enriched for genes that cluster with 
genes expressed in the germline and muscle/collagen, 
based on analysis of coexpression with genes known to 
function in certain tissues (Supplemental Fig. 2E ; Kim et al. 
2001). The subcellular localization of some of the proteins 
identified is known. Several reside in P granules, and many 
of the gene inactivations affect P granule structure, in- 
cluding csr-1 /Argonaute, rde-12/F58G11.2, adr-l/ADAR, 
puf-8 /pumilio, and npp-12/nucleoporin (Ariz et al. 2009; 
Updike and Strome 2009; Sheth et al. 2010). 

The ERI-1 /nuclease, ERI-6/7/helicase, and LIN-35/Rb 
are required for endogenous RNAi, negatively regulate or 
compete with exogenous RNAi, and are required for the 
somatic expression of many multicopy transgenes. To 
distinguish between gene inactivations that specifically 



silence the rol-6 ::lin-14 3' UTR transgene and gene in- 
activations that, like the eh genes, silence many multicopy 
transgenes, the 69 silencing gene inactivations were tested 
for silencing of a transgene bearing the rol-6(sul006) ::rol-6 
3' UTR and of five other GFP fusion transgenes with 
different promoters driving expression in diverse tissues 
and under the control of different 3' UTRs (Fig. 2B ; Sup- 
plemental Table 1). Most gene inactivations (60 out of 69 
tested) silence multiple transgenes, suggesting that they 
identify new eh genes, not negative regulators of miRNA 
function that could have emerged from this screen that 
used the rol-6 ::lin-14 3' UTR reporter gene (Fig. 2C; 
Supplemental Table 2). The nine gene inactivations that 
specifically silence the rol-6 ::lin-14 3' UTR reporter gene 
could identify negative regulators of miRNA function or 
factors that act in pathways that compete with an miRNA 
pathway for common factors, analogous to how the eh 
mutations enhance RNAi by debilitating endogenous 
silencing pathways. The few gene inactivations that also 
silence the rol-6(sul006) ::rol-6 3' UTR but not other 
transgenes could be roi-6-specific, such as the carboxy- 
peptidase gene suro-1 (Fig. 2C; Kim et al. 2011). 

To identify factors that, like the helicase gene eri-6/7, 
negatively regulate exogenous RNAi, we analyzed all 69 
gene inactivations for enhanced exogenous RNAi (Eri) 
pheno types. In two types of RNAi assays, wild- type 
worms were exposed to dsRNA of one of the 69 newly 
identified genes from the screen and dsRNA of several 
tester genes,- i.e., genes for which RNAi in wild type 
induces a weaker phenotype than the known loss-of- 
function phenotype (Supplemental Fig. 3A). If the candi- 
date gene inactivation causes enhanced RNAi response, 
RNAi of the tester gene at the same time would cause 
a more severe phenotype than the control. We identified 
22 eri genes by RNAi (Fig. 2B,D ; Supplemental Table 3). 
These include known eri genes such as the helicase gene 
eri-6/7 (Fischer et al. 2008) and the Argonaute gene ergo-1 
(Yigit et al. 2006) but also many genes that have not been 
previously implicated in RNAi. In parallel to the double 
RNAi experiments, 27 available viable mutants corre- 
sponding to a subset of the 69 gene inactivations that 
emerged from the genome screen were tested for en- 
hanced RNAi (Eri) phenotypes. Fourteen mutants, about 
half of the strains tested, display enhanced RNAi pheno- 
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A Conservation and essential functions of genes identified 



conservet 
in human 
(17 genes) 

conserved in 
human (52 genes) 




ERGO-1 endogenous 
RNAi pathway 



iin-54 DRM/dREAM complex 



: : 
usp-5 

aps-2 IAP2 complex 



C Silencing of multiple transgenes 
56 
silence 
multiple transgenes 



4 ro/-6-specific 
(incl. suro-1) 



rol-6::lin-14 
-specific 



D Defects in exogenous RNAi 



Enhanced 
RNAi (14) 

Enhanced RNAi 
(soma only) (8) 



Not 
determined 




RNAi-del 



ifective (6) 



x not Enhanced RNAi 
or RNAi-defective (2) 

siR-1/linc-22 silencing defects 



I Gene inactivation 


Description 


siR-1 target e 
20 °C 


xpression 
25 °C 


ergo-1 


Argonaute 


desilenced 


desilenced 


rde-4 


DsRBD 


desilenced 


desilenced 


eri-12/B0001.6 




desilenced 


desilenced 


drh-3 


helicase 


0 


desilenced 


mel-32 


SHMT 


desilenced 


desilenced 


mes-2 


SET EZH2 


0 


desilenced 


htz-1 


histone H2A.Z 


0 


desilenced 


puf-8 


Pumilio 




desilenced 


mes-3 


PRC2 complex 


0 


desilenced 


tcer-1 


TF 




desilenced 


rde-12 


helicase 


0 


desilenced 


F57C9.7 




0 


desilenced 


T03F6.4 




0 


desilenced 


Y39G10AR.7 




0 


desilenced 



Exogenous and endogenous RNAi defects 




Figure 2. A whole-genome screen to identify 
small RNA pathway factors. (A) Most of the 
69 genes identified are conserved in humans. 
Most of the identified genes are essential for 
viability. (B) Hierarchical clustering of the 
genes identified in the screen based on phe- 
notypes in RNA silencing assays. The clus- 
ters described in Figure 3A are indicated in 
color: Cluster 1 is burgundy, clusters 1 and 2 
are red, cluster 2 is orange, cluster 3 is dark 
blue, and cluster 4 is light blue. (C) Most gene 
inactivations silence multiple transgenes; 
some gene inactivations are specific for the 
rol-6 gene or the lin-14 3' UTR, representing 
possible rol-6 interactors or enhancers of 
miRNA activity, respectively. (D) Many en- 
hanced RNAi (eii) gene inactivations were 
identified, as well as several RNAi-defective 
gene inactivations (ide). [E] Genes required for 
the silencing of a single-copy siR-1 target 
transgene. Orange indicates moderate desi- 
lencing when the gene is inactivated, and red 
indicates strong desilencing. [F). Overlap of 
genes acting in exogenous RNAi and endoge- 
nous siR-1 silencing. 



types (Fig. 2B,D ; Supplemental Table 3), many of which 
confirmed the results of the double RNAi assays. 

Five of the genes identified are known to be required 
for RNAi: rde-4/ dsRBD, csi-1 /Argonaute, drii-3/heiicase, 
ekl-1 /Tudor, and htz-1 /H2A.Z. Interestingly, rde-4 mu- 
tants are mostly but not completely RNAi-defective but 
display enhanced RNAi phenotypes when inactivated by 
RNAi. Similarly, whereas a transgene crossed into a strain 
with an rde-4 mutation is initially silenced (e.g., the 
mgls30 transgene), after multiple generations, the trans- 
gene is no longer silenced. These data suggest that rde-4 
has multiple distinct functions, possibly in the soma 
and germline. We tested available viable mutants for 
RNAi-defective (Rde) phenotypes and assayed all 69 
genes for Rde phenotypes by RNAi. We identified one 
new factor required for RNAi, a DEAD-box helicase 
protein homologous to Drosophila Vasa and Belle and 
human DDX3 and DDX4 that localizes to germline 



ribonucleoprotein P granules (Sheth et al. 2010). We 
named this gene rde-12. 

From the suite of secondary screens, 28 out of 69 gene 
inactivations show altered responses to exogenous RNAi: 
either enhanced silencing in response to dsRNA (22 eh 
genes) or a lack of silencing in response to dsRNA (six rde 
genes). 

Endogenous RNAi defects caused by the newly 
identified gene inactivations 

Several known endogenous RNAi pathway factors were 
identified in the full genome screen for enhanced trans- 
gene silencing, including the genes encoding the heli- 
cases ERI-6/7 and DRH-3, the Argonautes ERGO-1 and 
CSR-1, and the Tudor domain protein EKL-1. Many of the 
currently known negative regulators of exogenous RNAi 
(eh genes) are required in the ERGO-1 endogenous RNAi 
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pathway that targets newly duplicated genes (Fischer 
et al. 2011). This pathway may compete for limiting 
factors with the exogenous RNAi pathway and transgene 
silencing pathway,- thus, loss of this endogenous silencing 
pathway may free factors to in turn enhance exogenous 
RNAi and transgene silencing. To test specifically for 
a role in endogenous RNAi, we assayed each gene in- 
activation for defects in endogenous RNAi silencing. We 
used a sensor strain that assays for the silencing activity 
of one abundant siRNA, siR-1, depending on the ERGO-1 
pathway and downstream factors that act in multiple 
endogenous RNAi pathways (Montgomery et al. 2012). 
siR-1 targets the noncoding RNA linc-22. The sensor 
strain harbors a gfp transgene expressed under the control 
of the strong ubl-1 promoter and the ubl-1 3' UTR with 
an insertion of a sequence that is the reverse complement 
of siR-1. The silencing of this transgene is dependent 
on siR-1, and an otherwise identical transgene without 
the siR-1 target sequence is not silenced (Montgomery 
et al. 2012). The siR-1 endogenous siRNA is abundant in 
C. elegans grown at 20°C but is 10-fold reduced in animals 
reared at 25°C (Supplemental Fig. 3B), allowing for more 
sensitive screening for factors required for endogenous 
siRNA-mediated silencing. We identified several new gene 
inactivations that disrupt silencing of this transgenic siR-1 
target but do not affect expression of the transgene that 
lacks the siR-1 target sequence (Fig. 2B,E), including the 
newly identified RNAi factors eri-12 (see below) and 
rde- 12 /DTLAD-box helicase. Many of the genes required 
for endogenous RNAi also affect exogenous RNAi (Fig. 2F). 

The identified RNAi pathway genes cluster 
into transcriptionally coregulated groups 

To study transcriptional coregulation of this set of 69 
genes, we extracted for each identified gene the 100 genes 
most similar in expression profile using SPELL (Serial 
Pattern of Expression Levels Locator), a search engine 
that uses —150 RNA expression data sets to find genes 
with similar expression profiles (Hibbs et al. 2007). 
Noting extensive overlap of coexpressed genes between 
the hits, the 69 genes were clustered using hierarchical 
clustering based on the overlap between the coexpressed 
genes (Fig. 3A). Most genes fall into four clusters, most 
notably a cluster (cluster 1) that contains nine genes, 
including the Argonaut e/piwi genes csr-1 and prg-1, the 
pumilio gene puf-8, and the succinate dehydrogenase 
gene sdha-2. Within this cluster, about half of the genes 
that are coexpressed with each of the hits from the screen 
are identical. The genes mostly frequently (up to 1 7 times) 
found coexpressed with the genes identified in the screen 
(Supplemental Table 4) include two Argonautes, transcrip- 
tion factors, chromatin factors, splicing factors, nucleo- 
porins, and RNA-binding proteins. For example, the attf-2 
AT hook transcription factor is coregulated with cluster 2 
and 3 genes and the Zn finger transcription factor somi-1 
that regulates a nuclear response to the miRNA mir-84 
(Hayes et al. 2011) with cluster 4 genes. 

To understand how these genes are coregulated, we 
examined published microarray data and ChIP data and 



looked for promoter motifs. We found an overlap between 
genes regulated by the interacting proteins MEP-l/Zn 
finger and LET-418/Mi-2 and cluster 2 genes,- -50% of 
cluster 2 genes are enriched in ChIP experiments for 
MEP-1. The expression data sets that contribute most to 
the similarity in expression profiles for each of the 
clusters include germline expression profiles and Mi-2/ 
NuRD complex mutants (cluster 1), LIN-35/Rb mutants 
(cluster 3), and DRM complex mutants and RNAi path- 
way mutants (cluster 4). Motif searching confirmed the 
connection between enhanced silencing and the Rb 
pathway (Wu et al. 2012). We used each of these four 
clusters to extract the 75 genes most similar in expression 
profile to these gene clusters. From these genes, we 
searched the sequences 600 base pairs (bp) upstream for 
motifs using MEME and identified several upstream 
motifs (Fig. 3B ; Supplemental Fig. 2F). Motif 1 was recently 
identified in the promoters of germline and oogenesis 
genes (Linhart et al. 2012) and ribosomal protein genes 
(Sleumer et al. 2012) as part of a bipartite motif, although 
the second element is not always present in cluster 1 
promoters. This second element of motif 1 corresponds to 
a germline E2F-binding site (Fig. 3B ; Kudron et al. 2013). 
Cluster 1 motif 2 is often found in close proximity to motif 
1 (Supplemental Fig. 2F). We also identified an E2F motif in 
cluster 2 genes (Fig. 3B ; Chittenden et al. 1991). E2F con- 
stitutes the DNA-binding protein of the retinoblastoma/ 
dREAM complex, many subunits of which emerged from 
the gene inactivation screen. Taken together, these data 
suggest that many of the genes identified are regulated by 
Mi-2/NuRD-like and DRM/dREAM complexes. 

eri-12 and ntl-9/CNOT9 are required for endogenous 
siRNAs 

Multiple endogenous RNAi pathways are active in 
C. elegans. To place the newly identified RNAi factors 
into specific endogenous RNAi pathways that are char- 
acterized by specific Argonaute proteins, we subjected 
the 69 gene inactivations to a series of genetic and cell 
biological assays (Fig. 2B ; Supplemental Table 2). More 
than half (40 of 69) of the gene inactivations show 
phenotypes indicative of defective endogenous RNAi 
pathways. Several of these were selected for further 
functional analysis, including the conserved genes arl-8/ 
ARL8 and ntl-9/CNOT9 as well as eri-12 /B0001. 6, a gene 
that clusters with ERGO-1 pathway genes in our pheno- 
typic analysis (Figs. 2B ; 4A). 

The gene B0001.6 encodes a protein of unknown func- 
tion with orthologs in nematodes. Consistent with a func- 
tion in RNAi, the B0001.6 protein was identified as an 
interactor of the RNAi pathway gene RDE-10 by tandem 
mass spectrometry analysis of proteins that immuno- 
precipitate with RDE-10 (Zhang et al. 2012). We obtained 
two B0001.6 mutants: a transposon insertion that disrupts 
a splice site required for splicing of the 5' UTR to the first 
exon (Supplemental Fig. 3C) and a missense mutation, 
gk214100 (G206E). Both of these mutants recapitulate 
the enhanced RNAi phenotype observed in double RNAi 
experiments (Fig. 4B) ; hence, we named this gene eri-12. 
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K CA ACCC--19ntor35nt-TCCCmAA 
T C CCTTTAA 



Cluster 1 Motif 1 

(4.3e-122; in 46 sites in 41/72 promoters within 600 bp of ATG) 

Germline, oogenesis and ribosomal protein genes (Linhart et al 2011, Sleumer et al 2012); 
prg-1; pgl-1; sdha-2 (17/72 cluster 1 genes) 

Germline E2F motif (Kudron et al 2013) 



j\j\ A Cluster 2 Motif 1 
9§§i (8.4e-034; 47 sites in 29/49 promoters within 500 bp of ATG); isw-1, htz-1, lin-9, lin-35 

C C AAA Consensus Mammalian E2F binding site (Chittenden et al 1991) 



Figure 3. The 69 genes identified in the screen cluster based on overlap in coexpressed genes. (A) Hierarchical clustering of genes 
coexpressed with each of the 69 genes identified in the screen. Indicated in color is the percentage of overlap between the 100 genes 
most similar in expression profile for each pair of genes identified in the screen. The clusters of genes that share many coexpressed 
genes are indicated. [B) Two motifs identified in cluster 1 and cluster 2 genes and in the genes most similar in expression profile to these 
genes. 



Loss of the gene eh- 12 causes a strong siR-1 sensor 
desilencing phenotype similar to the phenotype observed 
after inactivation of the known ERGO-1 pathway com- 
ponents (Fig. 2E). To determine whether eh- 12 is required 
for silencing of target mRNAs of the ERGO-1 pathway, 
we assayed for levels of target genes of endogenous 
siRNAs in the eh- 12 mutant. Using qRT-PCR, we found 
that the eh- 12 mutant fails to completely silence ERGO-1 
pathway targets (Fig. 4C). Using qRT-PCR, we also found 
that the eh- 12 mutant has reduced levels of endogenous 



siRNAs that depend on ERGO-1 (Fig. 4D). These results 
show that eh-12 acts in the ERGO-1 endogenous RNAi 
pathway and is required for endogenous siRNA biogenesis. 
Downstream from siRNA production, the Argonaute pro- 
tein NRDE-3 associates with ERGO-1 pathway siRNAs 
and translocates to the nucleus to mediate nuclear RNAi. 
A loss of ERGO-1 pathway siRNAs causes a shift in 
NRDE-3 : : GFP fluorescence from the nucleus to the cyto- 
plasm (Guang et al. 2008). Loss of eh-12 results in a re- 
duction in nuclear NRDE-3 : : GFP signal (Fig. 4E), indicat- 
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ing a loss of siRNAs binding to NRDE-3 to direct H3K9 
trimethylation. 

Several other newly identified genes appear to act in 
the ERGO-l-ERI-6/7 endogenous RNAi pathway, includ- 
ing C43H6.4, a gene conserved in nematodes encoding 
a protein with an Armadillo-type fold. Like eri-12, loss of 
C43H6.4 results in a reduction of ERGO-1 pathway 
silencing (Fig. 4C) caused by a loss of ERGO-l-ERI-6/7- 
dependent siRNAs (Fig. 4D). eri-12, C43H6.4, and nrde-3 
share expression profiles with an overlapping set of genes, 
indicative of a potential interaction between these genes 
(Fig. 4F). Similar to eri-12 and C43H6.4, inactivation of 
the gene ntl-9, one of several CCR4/NOT complex genes 
identified (Supplemental Table 8), results in defects in the 
ERGO-l-ERI-6/7 endogenous RNAi pathway (Fig. 4B-D). 
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Figure 4. Characterization of new factors in the ERGO-l-ERI- 
6/7 endogenous RNAi pathway. [A) Selected genes and their 
reduction-of-function phenotypes in assays for different RNAi 
pathways. Indicated in yellow are genes described in more detail 
here and in Figure 5. [B] eri-12 and ntl-9 inactivation enhances 
exogenous RNAi efficacy in RNAi assays targeting the genes 
hmr-1, dpy-13, unc-73, and lin-1. (C) eri-12, ntl-9, and C43H9.6 
affect expression of endogenous siRNA targets (F55C9.3 and 
F39E9.7) as determined by qRT-PCR. (D) eri-12, ntl-9, and 
C43H9.6 affect endogenous siRNA (22G siR-1 and 26G-01 
[C40A11.10]) levels as measured by qRT-PCR. (£) Loss of eri- 
12 causes a reduction of nuclear localized NRDE-3/ Argonaute. 
[F) Overlap between genes transcriptionally coregulated with 
eri-12, C43H6.4, and nrde-3. 



The ARF-like GTPase ARL-8 acts in the multiple 
endogenous RNAi pathways 

Using genetic interaction analyses, we assayed the gene 
inactivations for a role in endogenous RNAi pathways 
other than the ERGO-1 pathway (Fig. 4A). Genes acting in 
the CSR-1 endogenous RNAi pathway are known to 
interact with Ras pathway genes. We tested for enhance- 
ment of the ksr-l(n2526) mutant and found several gene 
inactivations (including arl-8) that enhance the pene- 
trance of LI lethality of this mutant (Fig. 4A). The arl-8 
gene scored positive as a negative regulator of exogenous 
RNAi in both double RNAi assays and tests of muta- 
tions in the gene, tm.2412 and tm.2388. arl-8 encodes 
a conserved small GTPase that functions in endosome-to- 
lysosome transport. In the genome-wide screen, we iden- 
tified several other endocytosis genes (including dpy-23, 
apa-2, and aps-2) that function in clathrin-dependent 
endocytosis. We rescreened a small set of endocytosis 
genes for transgene silencing using a protocol that uses 
a shorter exposure to dsRNA to allow for the identifica- 
tion of lethal genes. From this screen, we identified rme-4 
and rab-7 (Fig. 5 A), encoding another GTPase acting in 
endosome-to-lysosome trafficking. 

Because most of the endocytosis genes are essential, we 
focused on arl-8, for which viable reduction-of-function 
alleles are available (Klassen et al. 2010). The enhanced 
RNAi phenotype suggested that arl-8 functions in RNAi. 
Using deep-sequencing analysis, we determined the levels 
of the siR-1 siRNA, which is dependent on most known 
RNAi pathway genes, and found that this siRNA is re- 
duced by 50% (Fig. 5B). In an analysis of all known classes 
of C. elegans siRNAs, we found that antisense siRNAs 
targeting ERGO-l-ERI-6/7 target genes and ALG-3/4 tar- 
get genes are reduced by 25%-40% in the nonnull arl-8 
mutant (Fig. 5C; Supplemental Fig. 4A,B,D). Two types of 
siRNAs act in the ERGO-1 and ALG-3/4 pathways: 26G 
primary siRNAs and 22G secondary siRNAs. arl-8 affects 
26G siRNA levels: More than half of the abundant 
(>1 RPM) 26G siRNAs are more than threefold reduced 
as compared with <5% of non-26G siRNAs in the mutant 
(Fig. 5D). ERGO-1 target genes that produce higher levels 
of siRNAs are more affected by reduced ARL-8 (Supple- 
mental Fig. 4B). Whereas loss of ERI-6/7 results in a general 
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Figure 5. ARL-8 acts in endogenous RNAi. (A) Inactivation of endocytosis genes increases transgene silencing. Indicated in orange are 
gene inactivations that enhance somatic transgene silencing, "n/a" indicates that the phenotype could not be scored because of lethality. 
chc-1 /clathrin (RNAi) is lethal, but a temperature-sensitive chc-1 mutant displays Eri phenotypes. [B] Reduced levels of siR-1 in an arl-8 
mutant as determined by deep sequencing. (C) arl-8 mutation affects endogenous siRNAs in the ERGO-1 pathway and the ALG-3/4 
pathway. [D] arl-8 is required for most 26G siRNAs but not for 22G siRNAs. (E) arl-8 is required for siRNAs of ALG-3/4 target genes that 
produce mostly 26G siRNAs. (F) A model of how the different silencing pathways and factors identified affect transgene silencing. 



moderate reduction of ALG-3/4 target siRNA levels (Fig. 
5C ; Supplemental Fig. 4C), a subset of ALG-3/4 target 
genes is strongly affected in arl-8 mutants (Supplemental 
Fig. 4D). We found that ALG-3/4 target genes fall into two 
classes: One class of genes is targeted by a single 26G 
siRNA either close to the start codon or targeting the 3' 
UTR, and other ALG-3/4 target genes are targeted by 
multiple siRNAs distributed over the gene, including 
secondary 22G siRNAs. The ALG-3/4 target genes most 
depleted for siRNAs in arl-8 are genes that mostly produce 



26G siRNAs (Fig. 5E). We found that siRNAs corresponding 
to a subset of ALG-3/4 target genes are enriched in CSR-1 
immunoprecipitation experiments (Fig. 1H; Claycomb et al. 
2009). Thus, CSR-1 may act in parallel with or downstream 
from ALG-3/4. arl-8 affects siRNAs of ALG-3/4 targets that 
are not routed to the CSR-1 Argonaute (Supplemental Fig. 4E). 

arl-8 encodes a small GTPase that functions in endo- 
some-to-lysosome transport. In plants, ARL8 is also an 
essential host factor for tomato mosaic virus RNA 
replication that occurs at membranes. It interacts with 
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the viral superfamily I helicase and RdRP protein. Simi- 
larly, ARL-8 may interact with the ERI-6/7 superfamily I 
helicase (with which it shares a set of coexpressed genes) 
(Fig. 3A ; Supplemental Fig. 4F) and the RRF-3 RNA- 
dependent polymerase to promote 26G siRNA biogenesis 
or stability. 

Discussion 

Transgenes in many organisms are silenced and produce 
siRNAs that mediate RNAi-based silencing. Like trans- 
posons and viruses, upon cell entry, transgenes may be 
recognized as foreign and targeted by the RNAi pathway. 
However, the silencing continues every generation after 
the initial injection of transgenes, as if a "foreign" signal is 
continually regenerated. Our deep RNA sequencing anal- 
ysis revealed that only particular segments of the trans- 
gene are surveilled in this way, with siRNAs particularly 
concentrated in the lin-14 3' UTR of the fusion gene. One 
interpretation of the 3' UTR location of many of the 
transgene siRNAs is that the ribosome "clears" the mRNA 
of siRNAs in the coding region. Another interpretation is 
that endogenous miRNAs that target these 3' UTRs may 
attract the synthesis of secondary 22G siRNAs in a manner 
similar to how primary siRNAs trigger the production of 
secondary siRNAs, although the normal silencing that we 
observed in the absence of miRNA function did not sup- 
port this hypothesis. In fact, we found that loss of the 
miRNA factors ALG-l/Argonaute, AIN-1/GW182, and 
VIG-1 causes increased transgene silencing. It is possible 
that the miRNA machinery competes with the siRNA 
machinery for the lin-14 3' UTR or another limiting com- 
mon factor. Alternatively, ALG-1 may bind 3' UTR trans- 
gene siRNAs generated by Dicer and silence the transgene 
post-transcriptionally. VIG-1 was found in immunoprecip- 
itation experiments for the RNAi factor RDE-10, suggest- 
ing that miRNA and RNAi factors may interact. 

The production of siRNAs from the foreign transgene 
depends on the antiviral silencing factors RDE-1 /Argo- 
naute and DRH-1 /helicase and paradoxically on ERI-6/7/ 
helicase, a gene inactivation that causes increased trans- 
gene silencing. Silencing also depends on the nuclear 
factors NRDE-3/Argonaute, LIN-61, and SET-25/HMT. 
Silencing of the transgene in the ERI-6/7/helicase mutant 
is accompanied by changes in transgene chromatin mod- 
ifications: an increase in the trimethyl-H3K9 mark and 
a decrease in the trimethyl-H3K4 mark. This suggests 
that transgene siRNAs are bound by NRDE-3 to direct the 
trimethylation of histone H3 at K9 to cause transcrip- 
tional silencing. The reduction in transgene siRNAs in 
the eri-6/7 mutant may be a result of transcriptional gene 
silencing so that there is a resulting reduction in the 
mRNA template for siRNA generation. Alternatively, the 
ERI-6/7 helicase may be directly involved in the bio- 
genesis or stability of a subset of the transgene siRNAs. 
Finally, transgene siRNAs could be destabilized upon loss 
of the ERGO-l-ERI-6/7 endogenous RNAi pathway. 
Upon loss of the ERGO-l-ERI-6/7 endogenous RNAi 
pathway, more NRDE-3/ Argonaute may be available to 
bind transgene siRNAs that in wild type interact with 



other Argonautes. siRNAs of a subset of ERGO-l-ERI-6/7 
endogenous siRNA targets are bound by both NRDE-3 
and CSR-1, and we hypothesize that transgene siRNAs 
also are bound by multiple Argonaute proteins. Thus, 
CSR-1 may act in parallel with NRDE-3 or other WAGO 
proteins, downstream from RDE-1. Hence, in the absence 
of CSR-1, more transgene siRNAs may be bound by 
NRDE-3 to more effectively silence transgenes. 

Inactivation of dpy-30 and set-2 causes enhanced trans- 
gene silencing. Thus, the SET1/MLL complex is likely to 
deposit activating H3K4me3 marks on nonsilenced trans- 
genes, possibly directed by the anti-silencing Argonaute 
CSR-1. Other gene inactivations identified in our screen 
for enhanced transgene silencing may also act in pro- 
moting active transcription. NURF-1 is the ortholog of 
BPTF, a factor that binds H3K4me3 to recruit active 
transcription factors,- ELL-1 is an RNA polymerase elon- 
gation factor,- and TCER-1 and PUF-8 interact to post- 
transcriptionally maintain levels of certain germline 
mRNAs (Pushpa et al. 2013). CSR-1 was recently identi- 
fied as an interactor of FBF-1, a paralog of PUF-8, to 
attenuate translation elongation (Friend et al. 2012). 

Negative regulators of exogenous RNAi could be genes 
that act in competing endogenous small RNA pathways, 
as exemplified by the ERGO-l-ERI-6/7 pathway that 
silences duplicated genes (Fischer et al. 2011). Several 
newly identified genes negatively regulate exogenous 
RNAi and act in the ERGO-l-ERI-6/7 endogenous RNAi 
pathway, including C43H6.4, the CCR-4/NOT complex 
gene ntl-9/CNOT9, and eri-12. The ERI-12 protein forms 
a complex with the RNAi factor RDE-10. Besides canon- 
ical RNA-binding proteins such as the RDE-1 2/helicase, 
we identified several metabolic enzymes that have been 
identified as (candidate) mRNA-binding proteins in hu- 
man cells: SDHA-2/SDHA, MEL-32/SHMT2, and PDI-2/ 
P4HB (Castello et al. 2012). 

We found that several genes that function in endocy- 
tosis, specifically factors acting in clathrin-mediated 
endocytosis, are among the genes identified in the RNAi 
screen for enhanced transgene silencing. These are the 
gene inactivations that are the most pleiotropic in our 
screen, causing slow development and decreased brood 
sizes. However, they are specific in that many other gene 
inactivations that cause such pleiotropic developmental 
defects have no transgene silencing phenotype. The endo- 
somal pathway is required for RISC formation and turn- 
over, with RISC localizing to multivesicular bodies (MVBs) 
in addition to P bodies (Gibbings et al. 2009; Lee et al. 
2009). The endocytic pathway has also been implicated in 
the uptake of dsRNA in Drosophila (Saleh et al. 2006). We 
found that the <xr]-8/Arf-like GTPase is required for wild- 
type levels of 26G siRNAs in two endogenous RNAi 
pathways and may promote 26G siRNA biogenesis. 

Our screen for gene inactivations that enhance trans- 
gene silencing identified numerous additional small RNA 
pathway factors (Fig. 5F): factors acting in three endoge- 
nous RNAi pathways, factors enhancing exogenous 
RNAi, factors required for exogenous RNAi, and factors 
regulating expression of RNAi factors (Wu et al. 2012). 
Many of these genes are conserved in humans and are 
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essential for viability. A subset of the genes identified is 
expressed in the germ line and regulated by E2F, Mi-2/ 
NuRD complex factors, and DRM/dREAM complexes. 
Overlap of the genes coexpressed with each of the genes 
identified in the screen was used to identify clusters of 
genes that may act in distinct aspects of transgene 
silencing. Often, genes that encode proteins that act in 
one protein complex cluster together, such as CCR-4 and 
NTL-9; LIN-9, LIN-35, and LIN-54; DRH-3 and EKL-1; 
APA-2 and DPY-23; and ALG-1 and AIN-1. Also, genes 
that cluster together based on phenotype tend to cluster 
in the same coexpression clusters (Fig. 2B). Thus, we 
identified a large set of genes that act in distinct aspects of 
small RNA silencing. Our data also point to extensive 
interactions between small RNA pathways. 

Materials and methods 

C. elegans growth and handling; strains used 

All experiments were performed at 20°C unless stated otherwise. 
The genome- wide RNAi screen was done using strain GR1703 
[lim-6::gfp, col-10 r.lacZ ::lin-14, and rol-6(sul006) ::lin-14 
3' UTR]. For deep-sequencing analysis, see the complete list of 
strains in Supplemental Table 5. arl-8(tm2472) and arl-8(tm2388) 
are 3' UTR deletions, arl-8 mRNA levels are twofold reduced 
in tm2472 versus wild type. tm2388 has been shown to be a 
reduction-of -function allele (Klassen et al. 2010). 



Genome-wide screen for transgene silencing 

Embryos of strain GR1703 [mgls30) were isolated by treatments 
with bleach; five embryos were plated onto 12-well RNAi plates 
and seeded in duplicate with E. coli HT115-expressing dsRNA 
for each of the 18,371 genes represented in the Ahringer (Fraser 
et al. 2000) and Vidal (Rual et al. 2004) libraries. Plates were 
placed at 22°C. Fl progeny of these animals were scored at the 
L3/L4 larval stages for the percentage of animals that display the 
phenotype of a silenced transgene; i.e., a non-Rol phenotype. 
0: 0% non-Rol; 1: <20% non-Rol; 2: >20% non-Rol; 3: >50% non- 
Rol; and 4: >95 % non-Rol. Five-hundred-sixty-eight genes (3 % of 
C. elegans genes) with scores >2.5 were retested in four times. 
Ninety-three genes (16% of genes retested) retested in two out of 
four retests and three out of six total tests with an average score 
>2.5 (35% silenced) and were selected for further characteriza- 
tion. rol-6 interactors were excluded from further analysis. The 
remaining 69 RNAi clones were sequenced to confirm the 
identity of the gene targeted. The negative control (dsRNA not 
targeting a C. elegans gene) was scored 438 times and scored 
an average of 2% non-Rol with a maximum of 10% non-Rol. 
For two gene inactivations, the RNAi clone targeted two 
highly homologous duplicated genes: sams-3/sams-4 and adr-1/ 
D2005.1. For sams-3/4, we were unable to determine which gene 
inactivation was causal, since the DNA sequence identity is 100% 
for most of the gene, and therefore we indicate the hit as "sams-3/ 
4." For adr-1 and D2005.1, we analyzed mutants and observed 
silencing phenotypes for both and therefore included both. 

Secondary assays by RNAi 

All assays were done in at least duplicate. 

Genetic interaction with HEl006[rol-6(sul 006)) and GR1200(s<jt- 
l(scl3)) was tested by exposing these worms to the dsRNA of the 



93 strong positives. The genes that interacted with rol-6 and/or 
sqt-1 were excluded from further analysis, leaving 69 genes for 
further analysis. 

Transgene silencing assays were all done at 22°C ; these assays 
were done as described above for GR1703. Transgenes are 
described in Supplemental Table 1 . 

ksr-1 enhancement (Rocheleau et al. 2005) and lin-31 suppres- 
sion (Ding et al. 2005) were tested using feeding RNAi of the 69 
positives to ksr-l(n2526) and lin-31(nl053) mutants at 20°C and 
25°C, respectively, scoring the next generation for rod-like LI 
lethals and the presence of a single vulva or no vulva, respectively; 
ksr-2(RNAi) and ain-l(RNAi) were used as positive controls. 

To test for an enhanced RNAi response, two experimental 
setups were used in which the worms were exposed to RNAi by 
feeding of two genes: the Eri tester gene and each of the 69 
positives from the screen ("gene of interest 77 ). In sequential RNAi 
experiments, wild-type N2 worms were exposed first to the gene 
of interest from the LI stage until the next-generation LI stage,- 
subsequently, the LI worms were moved onto Eri tester plates 
for unc-73, hmr-1, and lin-1 and scored in the next generation as 
embryos [hmr-1) or young adults [unc-73 and lin-1). In simulta- 
neous RNAi-by-feeding experiments, worms were exposed to 
mixed E. coli RNAi-by-feeding lawns that contain a ratio of 4: 1 of 
gene of interest :Eri tester RNAi; LI larvae were placed on these 
lawns and scored in the next generation as embryos [hmr-1). 

In all enhanced RNAi assays and the siR-1 sensor desilencing 
assay, the phenotypes were compared with the positive control 
[eri-6/7) and the negative control (wild type) and reported relative 
to these controls. 

Mutant analyses of exogenous RNAi responses 

Enhanced RNAi (Eri) assays and assays for RNAi-defective (Rde) 
phenotypes were performed as described previously (Tabara et al. 
1999; Fischer et al. 2008) using feeding RNAi for lin-1, unc-73, 
dpy-13, unc-33, hmr-1, unc-22, andpos-L Animals were plated 
as starved LI larvae and scored as adults [hmr-1 and pos-1 for 
progeny) or in the next generation (other clones). All assays were 
done at least twice and in duplicate. In all enhanced RNAi 
assays, the phenotypes were compared with the positive control 
[eri-6/7) and the negative control (wild type) and reported relative 
to these controls. 

ChIP 

ChlP-qPCR was done as described (Guang et al. 2010) using ~72-h 
adults and anti-trimethyl-H3K9 (Lys9) antibody (Milllipore, cata- 
log no. 07-523) or anti-trimethyl-H3K4 (Lys9) (Abeam, catalog no. 
ab8580). The rol-6 primers are specific for the sul006 mutation. 

RNA isolation 

Large-scale total RNA was isolated by dounce homogenization 
in RNA-Bee (Tel-Test) followed by chloroform extraction and 
isopropanol precipitation. Small-scale total RNA isolations were 
done in Trizol using pestles in eppendorf tubes to homogenize 
the embryos/worms. 

qRT-PCR on endogenous RNAi target genes 

For qRT-PCR, total RNA from embryos and/or young adults was 
DNase-treated using the Turbo DNA-free kit (Applied Biosys- 
tems). cDNA was synthesized using RETROscript (Ambion) 
following the vendor's protocol. qPCR was done in triplicate 
with a CFX96 machine (Bio-Rad) using iQ SYBR Green Supermix 
(Bio-Rad). Relative mRNA levels were calculated based on the 
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2~AAct method using the gene Y45F10D.4 for normalization. Two 
or three biological replicates were done for each PCR. The error 
bars indicate the standard error between the biological repli- 
cates. Primer sequences are listed in Supplemental Table 6. 

TaqMan qRT-PCR to detect small RNAs 

TaqMan qRT-PCR and data analysis were carried out as de- 
scribed (Han et al. 2009). Total RNA was extracted from embryos 
and 52- to 56-h post-hatching L4/young adult worms. At least 
two biological replicates were done for each assay. The error bars 
indicate the standard error calculated for the biological repli- 
cates. The TaqMan assays are described in Supplemental Table 7. 

Small RNA sequencing 

Total RNA was isolated from embryos. Eighteen-nucleotide to 
28-nt small RNAs were size-selected, tobacco acid pyrophospha- 
tase-treated to remove 5 ' di- and triphosphate groups, ligated to 
5' and 3' adapters, and subjected to RT-PCR according to the pro- 
tocol by Gu et al. (2009). Small RNA amplicons were sequenced 
using an Illumina genome analyzer. 

Small RNA data analysis 

Data analysis was done as described (Shi et al. 2013). Briefly, 
small RNA sequences were parsed and mapped to the C. elegans 
genome (WormBase release WS230) and transgene sequences. 
Small RNA reads from each library were normalized to the total 
number of mapped reads. The numbers of small RNA reads for 
small RNA sequences mapping to multiple positions in the 
genome or transgene sequence were divided by the correspond- 
ing number of genomic loci. 

Bioinformatics 

FunCoup (Alexeyenko and Sonnhammer 2009) was used to 
analyze interactions between genes and gene products at a con- 
fidence cutoff of 0.10. Using SPELL (Hibbs et al. 2007), other 
positives or known small RNA pathway genes that fall within 
the first 100 genes with the most similar expression profiles were 
identified. Expression mountains (microarray "topology map" 
positions) were determined for each gene (Kim et al. 2001). 
P-values were calculated using hypergeometric distribution 
analysis. 

Accession numbers 

Small RNA Illumina deep-sequencing data are available at the Gene 
Expression Omnibus (accession nos. GSE52785 and GSE32366). 
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